%% Assignment 5.1

close all; clearvars; clc;

%% 5.1.a
% To saturate 100 cm3 dry soil, 40 gram of water has to be added to the
% dry soil. What is the porosity of the soil?

% Assuming no airentrapment
p   = 40/100;   % [-]

fprintf('The porosity is %f\n',p);

%% 5.1.b
% If one wants to saturate the soil of exc. 1a) to a depth of 25 cm, how 
% much rainwater (volume basis: cm3) has to fall per cm2 on the soil 
% surface? One can assume that no rainwater evaporates or be lost or 
% stored at the surface.

d   = 25; % [cm]

fprintf('One needs %f [cm] of water column to saturate 25 [cm] of soil.\n',...
    d*p);

%% 5.1.c
% If the amount of rainwater is added instantaneously. What is then the 
% thickness of the water layer at the surface?

fprintf(['At t=0, the entire water column floats on top of the soil. With'...
    'a thickness of %f [cm]. But I thinck the prof means the d given'... 
    'above, %f [cm]\n'],d*p,d);